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1. Introduction 

The goal of research in high-energy heavy-ion collisions is to study the properties 
of strongly interacting matter at extreme energy density, including the possible 
phase transition to a colour-deconfined state: the Quark Gluon Plasma (QGP). The 
matter produced in these collisions can be probed using hadrons produced in partonic 
processes with a large momentum transfer ('hard scatterings'). These processes take 
place early in the collision and are only sensitive to short distance scales. In the absence 
of nuclear effects the hard production yields in nucleus-nucleus collisions are therefore 
expected to scale as if the collision were an independent superposition of nucleon- 
nuclcon collisions. Measurements at SPS have shown that dilepton production in the 
Drcll-Yan process indeed follows this expectation £Q. 

Among the first measurements at RHIC was the measurement of light hadron 
production at high transverse momentum p^ in Au+Au collisions, which shows a 
suppression with respect to the scaled p+p results. Since these first observations, 
measurements in d+Au collisions have confirmed that the observed suppression is a 
final state effect. Recently there has also been an increased activity to verify that high- 
Pt light-hadron production in proton-proton collisions can be understood in terms of 
perturbative QCD (pQCD) calculations, as expected for hard processes. Some of the 
relevant results will be reviewed in the next section. 

While for light hadron production much of the groundwork has been done and 
analyses are clearly moving towards more advanced observables like identified hadron 
spectra and correlation measurements, first results on open charm production are 
becoming available. Like high-p^ light hadrons, open charm is expected to be 
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dominantly produced in hard processes and can therefore serve as a calibrated probe 
of the medium. Due to their large mass, however, charm quarks and hadrons are 
expected to be affected differently by the medium than high-p^ light hadrons 

In the second part of this paper an overview will be presented of the existing 
results on open charm production at RHIC. The present results are based on run-2 
and run-3 data. First results from the large statistics Au+Au data sample from run-4 
are to be expected soon. These will greatly improve the precision and py-coverage of 
the open charm measurements in Au+Au collisions at RHIC. 

2. High-pT light hadron production 

High-px hadron production at is the most readily accessible observable for hard 
processes at RHIC. At sufficiently high pt, all hadrons are expected to be produced 
in jet fragmentation. The non-perturbative dynamics of jet-fragmentation can be 
characterised by a universal fragmentation function D(z) which is parametrised using 
data from e + e~ collisions at different energies j3] 0|. With these fragmentation 
functions and the parton densities from deep inelastic scattering experiments, 
the expected cross sections for high-pr hadron production can be calculated in 
perturbative QCD (pQCD). 

2.1. Neutral pions and charged hadrons in p+p collisions 

Figure shows pT-spectra of 7r° measured by PHENIX jS] (left panel) and charged 
hadrons from STAR |S] and BRAHMS 7 (right panel) in p+p collisions at y/s = 
200 GeV. The data are compared to a next-to-leading order (NLO) pQCD calculation 
jHJ- The uncertainty in the theoretical calculation is estimated by varying the 
renormalisation and factorisation scales fin and fip to half and twice the nominal value 
of p,ji = ii f = Pt- These scale variations change the calculated cross-sections by about 
20% for px > 5 GeV. To illustrate the uncertainty in the fragmentation functions, 
the calculation was performed with two different sets of fragmentation functions, 
from Kniehl, Kramer and Potter (KKP) [2] and from Kretzer Both sets were 
independently determined from similar selections of e + e~ data using slightly different 
assumptions about relations between the fragmentation functions for different partons. 
This turns out to be the dominant source of uncertainty for the ir° spectrum: variations 
of the order of 50% are seen, mainly due to uncertainties in the gluon fragmentation 
function. Note, however, that these are partly normalisation uncertainties and do 
not change the shape of the spectra very much. The measurements have an overall 
normalisation uncertainty of about 10%, which is not indicated in the figures. The 
systematic offset between the STAR and BRAHMS results in figure ^ can probably 
be attributed to this normalisation uncertainty. 

Both for neutral pions and charged hadrons, data and theory agree over more 
than 5 orders of magnitude, which gives confidence that hadron production at high 
Pt (> 3 GeV) is indeed governed by hard point-like processes. 

2.2. Strange hadron production in p+p collisions 

The above comparisons can be extended to the strange hadrons Kg and A. While 
the kaon fragmentation function is relatively well-constrained by the data from e + e~ 
collisions, data on A production are scarce 




Figure 1. Measured py-distributions of n° (left) and charged hadrons (right) in 
p+p collision compared to NLO pQCD calculations (lines). The different lines 
indicate results using different paramctrisations of the fragmentation functions 
(see text) and choice of the factorisation and renormalisation scales. The 
lower panels show the difference between data and calculations, divided by the 
calculated curves. 




Figure 2. Comparison of Kg and A spectra in p+p collisions as measured by 
STAR to NLO pQCD calculations 1101 . The different lines indicate results for 
several settings of the factorisation and normalisation scales. 
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Figure 3. Nuclear modification factor Raa for tt° and charged hadron 
production in central and peripheral Au+Au collisions, as measured by PHENIX 
[12) and STAR 0. 



A first comparison of Kg and A spectra in *Js = 200 GcV p+p collisions to 
NLO calculations as presented at this conference is shown in figure [21 [H3- The 
agreement between the measured K s spectrum and the expectation from NLO pQCD 
is reasonable, although the shape of the calculated spectrum is slightly more concave 
than the measured one. This difference is mainly apparent at relatively low px (1-2 
GeV), where soft production processes may still contribute significantly. 

For the A on the other hand, the agreement between the data and the NLO 
calculations is far from satisfactory. This might be indicative of the breakdown 
of the massless formalism and the factorization ansatz for particles with mass that 
is significant compared to px j 1 11 - Before drawing this conclusion, however, the 
uncertainties in the A fragmentation functions should be better quantified. 



2.3. Suppression in Au+Au collisions 

To compare measured the particle spectra in Au+Au collisions to the expected N co u 
scaling from p+p, the nuclear modification factor 

dN/dp T \ Au+Au 

UAA N coll dN/dp T \ p+p {1) 

is generally used. In figure [3| these ratios are shown for peripheral and central Au+Au 
collisions at = 200 GeV, both for tt° from PHENIX (left panel) [E] and charged 
hadrons from STAR (right panel) p|. The suppression ratio Raa in peripheral 
collisions is close to unity, while for central collisions a suppression of up to a factor 
5 is observed. This suppression was one of the first indications of a strong final state 
modification of particle production in Au+Au collisions that is now generally ascribed 
to energy loss of the fragmenting parton in the hot and dense medium. 



3. Open charm production in d+Au and Au+Au collisions 

While light hadron production is only expected to be calculable in perturbative QCD 
at higher pt, the charm quark mass (m c « 1.35 GeV) is large enough to expect pQCD 
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Figure 4. (left) Collected px spectra for D°, and D*^ production in 

d+Au collisions at yfs = 200 GeV. The spectra have been divided by N co u ss 7 
to be comparable to calculations for p+p. The and D*^ spectra were 

scaled to match the D° data, (right) Energy dependence of the total cross- 
section for charm production 17 18 . The different symbols indicate different 
experiments and collision systems. All data were scaled to effective p+p cross 
sections. When available, the statistical and systematic error were added in 
quadrature. The bands in both panels show NLO pQCD calculations 1141 with 
two different assumptions for the charm mass and two choices of the factorisation 
and normalisation scales (see test), variations. 



calculations to be valid for all px- Final state effects on charm quarks and hadrons 
are expected to be smaller than for the light hadrons due to the large charm mass [2] ■ 

3.1. Charmed meson spectra in d+Au collisions 

None of the RHIC experiments currently has a vertex detector with sufficient resolution 
to reconstruct secondary vertices of charm decays. Even without secondary vertex 
reconstruction, STAR has been able to statistically reconstruct decays to charged 
hadrons in d+Au collisions at y/s = 200 GeV using an invariant mass method |13| . 
Different charmed mesons (D°, D , and D*^) can be reconstructed in different px 
ranges, thus providing a charm measurement up to = 11 GeV, as shown in the 
left hand panel of figure 01 The Z? ± (triangles) and D*^ (squares) spectra have been 
scaled to match the D° spectra (circles) . Note that the D° has been measured at low 
Pt and at high px using different decay modes and thus provides a normalisation for 
the whole px range. 

Also shown in figure |3| (left) are NLO pQCD calculations of charm quark spectra 
14 . To illustrate the sensitivity of those calculations to the charm quark mass m c 
and the choice of renormalisation scale, curves are drawn for m c — 1.2 GeV and 
m c = 1.5 GeV and with two choices for the renormalisation and factorisation scales: 
Mii = A*f = m T and fiR = fj,p = 2mx- Note that the shape of the spectra at low px 
is most sensitive to both the charm quark mass and the choice of scales. 

For a detailed comparison of the data to theory, the calculated charm quark 
spectrum should be convoluted with the charm fragmentation function. This would 
lead to a softening of the spectrum and a reduction of the yield, both of which may 
in principle depend on the meson species. Given the limited px range of the spectra 
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for the separate species and the relatively large uncertainties in their fragmentation 
functions, we have chosen to compare the shape of the combined meson spectrum 
directly to the charm quark spectra. For this purpose, the calculated spectra were 
scaled up by a factor of 4 to approximately match the data. The shapes of the 
calculated charm quark spectra and the measured meson spectra are surprisingly 
similar, leaving little room for softening due to fragmentation. 

All in all, it is far from clear that the present data can be matched with a 
NLO pQCD calculation. Before drawing any conclusions about charm production 
in elementary collisions at RHIC, though, we should wait until the present data are 
finalised. Although it is expected that N co u scaling is valid for charm production in 
d+Au, it would be good to confirm this by similar measurements in p+p. There are 
also some open questions for theory. For example, a matched next-to-leading logarithm 
calculation is needed to describe beauty production at the Tevatron |15| . In addition, a 
new way of extracting fragmentation functions, by fitting the Mellin moments instead 
of a direct z-space fit is found to lead to an effectively harder fragmentation function 
16 . Similar considerations may also affect calculations of charm production at RHIC. 

3.2. Total charm cross section 

The right-hand panel of figure 0] shows the measured energy dependence of the 
total charm quark cross section ^| ^], compared to NLO pQCD calculations 
|14|. Estimating the total charm quark cross section from experimental data 
involves substantial extrapolations to the unmeasured regions of momentum space 
and corrections to include unmeasured charmed hadron species. These corrections 
lead to sizeable systematic uncertainties on the data, as can be seen from the figure. 
The uncertainties on the NLO pQCD calculations due to higher order corrections and 
the choice of the charm mass are also significant. It is therefore preferable to directly 
compare data and calculations in the measured regions (see also |19|1. 

Note also that it seems that the charm cross section per nucleon-nucleon collision 
as measured by STAR in d+Au collisions from a combination of the electron 
measurements and the invariant mass method is somewhat higher than expected from 
the trend observed by other experiments and the pQCD calculations. The deviations 
are within the present uncertainties on the measurements. 

3.3. Centrality dependence of charm production in Au+Au collisions 

A first indication of the centrality dependence of charm production in Au+Au 
collisions can be taken from electron spectra. After subtraction of the contributions 
from light hadrons (mainly through photon conversions, but also from Dalitz decays of 
7T°, ?7, rj' , p, u> and 4>) the electrons from heavy flavour decays remain ('non-photonic' 
electrons). In the left hand panel of figure the electron spectra from heavy flavour 
decays as measured by PHENIX in centrality-selectcd Au+Au collisions 20 are shown. 
The lines show reference spectra obtained from a fit to the measured spectrum in 
p+p, scaled by the number of collisions. At each centrality, the spectra agree with the 
expected N co u scaling from p+p, albeit within large errors. 

In the right-hand panel of figure |SJ the yields of non-photonic electrons with 
0.8 < pt < 4.0 GeV per nucleon-nucleon collision are shown as a function of centrality. 
There is no indication of a suppression as seen for light hadrons (see figure IHJl ■ One 
should keep in mind, however, that electrons with pt > 0.8 GeV have contributions 
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Figure 5. (left) Non-photonic electron spectra from Au+Au collisions, compared 
to the expectation from p+p. (right) Nuclear modification factors Raa for non- 
photonic electrons with px > 0.8 GeV as measured by PHENIX 1201 . 
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Figure 6. Elliptic flow V2 of non-photonic electrons in Au+Au collisions as 
measured by STAR and PHENIX |2l]. 



from semi-leptonic charm decays at all pr- The presented results are therefore not 
very sensitive to a possible suppression of charm production at moderate or high px 
(> 2 GeV). 

3.4- Charm flow 

A measurement of the elliptic flow t>2 of charmed mesons is an independent way of 
assessing the sensitivity of charmed mesons to final state interactions. Here again, we 
have to rely on measurements of decay electrons for the time being. In figure [B] the 
elliptic flow of non-photonic electrons is shown . Both STAR and PHENIX observe 
non-zero electron flow, which is a strong indication that charmed mesons flow. This 
is an intriguing possibility, because it would show decisively that charm production 
is sensitive to the dense hadronic or even partonic environment in the collision. At 
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the moment the statistical and systematic errors are still large, precluding a precise 
quantitative extraction of flow values. The situation is expected to dramatically 
improve with the larger Au+Au data samples which were recorded this year. 

4. Summary and outlook 

A comparison of neutral pion and charged hadron pT-spectra measured in p+p 
collisions at y/s = 200 GeV at RHIC to NLO pQCD calculations shows that high-pr 
light hadron production is well described by perturbative QCD, albeit within relatively 
large uncertainties, mainly from the fragmentation functions. This gives confidence 
that high-pr particle production is governed by hard, point-like processes, for which 
the cross section in Au+Au collisions is expected to scale with the number of nucleon- 
nucleon collisions. A suppression of light hadrons by approximately a factor of 5 
compared to the expected scaling is observed in central Au+Au collisions, due to final 
state interactions of the fragmenting quarks and/or the produced hadrons. 

For strange hadrons, Kg and A, the agreement between data and NLO pQCD 
calculations is not as good, or even unsatisfactory (for A). In those cases, however, 
the data do not extend to very high px and the fragmentation functions are not as 
well known as for the light hadrons. This needs more investigation before conclusions 
can be drawn about the applicability of pQCD. 

The same is true for open charm production in d+Au collisions, where shape 
of the measured D meson spectra is similar to the calculated charm quark spectra, 
leaving little room for softening due to fragmentation. 

First results on electron production from PHENIX indicate that there is no or 
very little suppression of charm production in Au+Au collisions. Measurements of 
electron flow, on the other hand, indicate significant flow of the charmed mesons, 
which can only be due to significant final state interactions. 

In the near future, a measurement of the p-r-dependence of nuclear modification 
factors for non-photonic electrons or maybe even open charm can be expected from 
the large statistics Au+Au data samples collected in run-4 at RHIC. This, together 
with a more accurate measurement of charm flow, will map out the interactions of 
charm quarks and hadrons with the medium and, through comparison with the light 
hadron results, may eventually shed more light on the nature of these interactions. 
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